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I Abstract 

7—i . A random polymer model is a one-dimensional Jacobi matrix randomly composed of 

^ I two finite building blocks. If the two associated transfer matrices commute, the corre- 

^ ' sponding energy is called critical. Such critical energies appear in physical models, an 

■ example being the widely studied random dimer model. It is proven that the Lyapunov 

I exponent vanishes quadratically at a generic critical energy and that the density of states is 

' positive there. Large deviation estimates around these asymptotics allow to prove optimal 

' lower bounds on quantum transport, showing that it is almost surely overdiffusive even 

, though the models are known to have pure-point spectrum with exponentially localized 

eigenstates for almost every configuration of the polymers. Furthermore, the level spacing 

• is shown to be regular at the critical energy. 

> ■ 
■ 

fi: 

> : 

X 1 Introduction 

Until quite recently, common wisdom was that one- dimensional random Schrodinger operators 
are in a strong localization phase and that there is nothing else of any interest to be discovered. 
In 1990, Dunlap, Wu and Phillips [DWP^ studied the random dimer model. It is similar to the 
conventional Bernoulli-Anderson model with a Hamiltonian given by the sum of the discrete 
Laplacian and a random potential taking only two values, except that these potential values now 
always come in neighboring pairs (dimers). For suitable values of the parameters, this model has 
so-called critical energies at which the two transfer matrices across the dimers commute. This 
leads to a vanishing of the Lyapunov exponent and a divergence of the localization length at these 
energies. They then argued and showed numerically that the second moment of the position 
operator X on the lattice grows superdiffusively under the dynamics like {^\X'^(t)\'ip) ~ Ct^/"^ 
for a localized initial state if) and any typical dimer configuration. This is likely to be responsible 
for the high conductivities of certain organic polymer chains jPWi) and quasi-lD semiconductor 
superlattices |PTBj . 
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This work considers Jacobi matrices H^^ ramdomly composed of two building blocks modeling 
two different polymers. The associated polymer transfer matrices are supposed to commute at 
an energy lying within the spectra of the periodic operators containing only one of the polymers. 
It is proven that for almost every polymer configuration u and every a > there is a positive 
constant Ca such that 

^ (0|e*^*|X|''e-*^*|0) > C^T^-^-'^, (1) 

where |0) is the state localized at the origin. For q = 2, this rigorously confirms that the 
heuristics of [DWP. (discussed below) provide a correct lower bound on transport. To prove the 
corresponding upper bound for critical energies of order 1 in the sense of Definit ion |21 remains an 
open problem, but we believe the quantitative lower bound on the Lyapunov exponent (Theorem 
12)) to be a central ingredient. Moreover, it is shown that for every configuration the l.h.s. is 
greater or equal than C T'^^^ for some C > 0. Note that ((H) implies that the conductivity is 
infinite either at finite temperature or if the critical energy is at the Fermi level [SB] . 

The above results should be confronted with the fact that the spectrum of a random polymer 
model is almost surely pure-point with exponentially localized eigenfunctions. For the related 
Bernoulli- Anderson model, such spectral localization results were first proven in |CKMj . later 



on in |S VWj . More recently, the random dimer model and the continuous Bernoulli- Anderson 
model were treated in [BG, and USSlj respectively. These works also established dynamical 
localization on energy intervals not containing a discrete set of special energies which includes 
the above critical energies. While |DSSlj considers continuum models, its approach can be 



carried over to prove spectral localization and dynamical localization away from the set of 
special energies for the polymer models studied here, see |DSS2j for the more general case of an 
arbitrary number of building blocks of bounded length. 

The fact that spectral localization can in principle coexist with quantum transport (even 
almost ballistic; note that ballistic is impossible with pure point spectrum |Simj ) was demon- 
strated by an example in |RJLSj (see also |BTp . However, those examples were rather artificial 
and much research since then was devoted to the program of proving dynamical localization (i.e. 
boundedness in time of the left-hand side of (^) in models of physical interest with previously 
estabhshed spectral localization, by means of upgrading the proof of pure point spectrum (see 
|GKj and references therein, also |BJj ). The success of this program may have raised doubts 
as to the validity of the distinction between spectral and dynamical localization in physically 
relevant contexts. This paper demonstrates that the distinction should indeed be made as it 
shows that exponential localization and quantum transport coexist also in physical models. 

Let us sketch the heuristics of |DWPj leading to ([Q). It is known |Bovj and proven below that 
the Lyapunov exponent generically vanishes quadratically like 7(-E'c + e) = ce^ + 0{e^) in the 
vicinity of the critical energy Ec- The extension of the eigenstates in an e-neighborhood of Ec 
is given by their localization length equal to the inverse of the Lyapunov exponent. Therefore 
the portion of the initial wave packet lying energetically in this e-neighborhood spreads out 
ballistically up to time scales T ^ e~^. Because it will be shown that the density of states is 
positive at Ec, this portion of states is proportional to e. Consequently, the gth moment of 
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the position operator should grow hke T'^ multiphed by this factor e ~ T showing that (jT]) 
should hold with high probability. 

The main technical tools are adequate action-angle variables, also called modified Priifer 
variables in the mathematical literature. Adapting techniques from jPFj, one obtains pertur- 
bative expansions around the critical energy of both the density of states and the Lyapunov 
exponent, which prove positivity of the density of states and quadratic vanishing of the Lya- 
punov exponent near a generic critical energy. The proof of requires an additional large 
deviations analysis for the Lyapunov exponent. The methods of proof are calculatory, quanti- 
tative and optimal. For example, they allow to show how large the moments of the position 
operator have to be if the commutator of the transfer matrices is small, but does not vanish. 

Acknowledgements: This paper is a heavily revised version of a preprint |JSSj which con- 
tained a first, but less direct proof of the deterministic lower bound stated in Theorem^ below. 
The basic strategy (Lemma Hand Section IB]) of the present proof of the dynamical lower bound 
(Theorem |3]) given the boundedness of transfer matrices (Theorem Ej) was suggested by S. 
Tcheremchantsev. This technique, which will be published in full generality in jDT], is simpler 
than the Guarneri method |Guaj of proving lower bounds employed in |JSSj (and also applicable 
here) and allowed us to circumvent previous more intricate arguments. We greatly appreciate 
that S. Tcheremchantsev made his work available prior to publication. S. J. and H. S.-B. were 
supported by NSF grant DMS-0070755, H. S.-B. moreover by DFG grant SCHU 1358/1-1 and 
the SFB 288. G. S. was supported by NSF grant DMS-0070343. He would also like to acknowl- 
edge financial support of CNRS (France) and hospitality at Universite Paris 7, where part of 
this work was done. 

2 Model and main results 

Let £-1- = (t±(0), . . . , £-|-(L-|- — 1)) and v± = {v±{0), . . . ,v±{L± — 1)) be two pairs of finite 
sequences of real numbers, satisfying t±{l) > for all / = 0, ... , L± — 1. These numbers are the 
hopping and potential terms of two different polymers. A family of random Jacobi matrices is 
now constructed by random juxtaposition of these polymers. More precisely, to any sequence 
uj = {(jJi)i(zi of signs -|- and — one associates sequences C = (C('^))nez and v^j = {vi^{n))n£Z 
by means of tuj = {. . . ,tujo,t^-^, . . .) and v^^ = (. . . , -O^^y, -0^^, . . .). An exact definition of the 
underlying probability space (fi, P), which also requires to randomize the position of Vujo{0) and 
Co(0), is given in Section [4.11 The polymer Hamiltonian H^^ of the configuration u is then 
defined by 

{H^ij){n) = -t^{n + l)^(n + 1) + t;,(n)^(n) - t^(n)^(n - 1) , G ^^(Z) , (2) 

and {Hi^)^<=Q becomes a family of random operators if the signs are chosen with probabilities 
p+ and p_ = 1 — J9+ respectively. The polymer transfer matrices Tj' at energy E E are 
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introduced by 



^± - '^v±{L±-i)-E,i±{L±-i) ■ ■ ■T{,^{o)~E,t±{o) > where '^'"'^ ~ J I ) ' ^'^^ 
The transfer matrices over several polymers are then 

T^(k, m)=Tf -T^ -...-T^ , k>m , (4) 

and T^{k,m) = T^{m, k)~^ ii k < m, T^{m,m) = 1. The Lyapunov exponent at energy E, 
also called inverse localization length, is then almost surely defined by (some more details are 
given in Sections 13.41 and H^Tj) 

7(i?) = £rn ^^log(||rJ^(A:,0)||) , (5) 

where (c±) = p+c+ +]9_c_ for any complex numbers c±. Vanishing of the Lyapunov exponent 
is considered an indicator for possible delocalization. For a polymer chain, this happens in the 
following situation: 

Definition 1 An energy Ec E M is called critical for the random family {Hi^)^(=q of polymer 
Hamiltonians if the polymer transfer matrices T^" are elliptic {i.e. |Tr(Tj''')| < 2) or equal to 
±1 and commute 

[T5,Tfl = . (6) 

Remark 1 The definition does not allow the critical energy to be in a spectral gap or at the 
band edges of one of the periodic operators (constructed from (^4., -0+) and (t_, {)_) respectively) 
except for points of band touching (where the transfer matrix is ±1). 

Remark 2 The condition (jU)) contains 4 equations. Given a model, one can only vary the 
energy. Hence, in the space of polymer models existence of critical energies is a non-generic 
property. On the other hand, given an energy Ec, it is always possible to construct polymer 
models that have E^ as a critical energy. 

Examples If L± = 1, the model reduces to the Bernoulli- Anderson model and there are no 
critical energies. If = 2 and L_ = 1, an example can be constructed as follows: choose 
t{l) = 1 for all I E Z and -0+ = (0,0) and v_ = (A) and |A| < 2, then Ec = is the critical 
energy. The most prominent jDWPl \Bov\ IBGj example is the random dimer model for which 
L+ = L_ = 2 and v+{0) = v+{l) = X and V-{0) = v^l) = -A (A G M), and t{l) = 1 for all 
/ G Z. This model has two critical energies Ec = X and Ec = —X as long as A < 1. It was 
previously (non-rigorously) known that 7(-Ec + e) = C^(e^) |Bovj for the random dimer model. 

The definition of the critical energy assures that there exists a real invertible matrix M 
transforming T^" and T^" simultaneously into rotations by angles ?7_ and rj^ respectively: 

V sm(r7±j cos(?7j-1 / 



4 



Hence 7(-E'c) = 0. Because are polynomials of degree L± in E, one can expand T^"'^'^ into 
powers of e. Since MTf^M'^ are rotations, this implies that ||MTJ''''''^M~-^|| < 1 + c|e| for 
|e| < Co and one deduces the following: 

Proposition 1 For eg > there exists a constant C < oo such that for all \e\ < eg and 
m,k E Z, 

\\T^^+'{k,m)\\ < C'e^l^ll'^-™! . (8) 
In particular, |7(-Ec + e)| ^ C" |e| for C > 0. 

Note that the bound in (jH)) does not depend on the configuration. To study a possible 
spreading of wave packets due to the divergence of the localization length, one best considers 
the moments of the associated probability distribution, notably the time-averaged moments of 
the position operator X on £^(Z): 

M^,,{T) = -e-^ (Ole^^-^lXre-^^-^lO) , g>0, (9) 
Jo ^ 

The exponential time average may be replaced by a Cesaro mean without changing the asymp- 
totics (e.g. )GSBj ). Proposition [H will lead more or less directly to the following deterministic 
lower bound on transport. 

Theorem 1 There exists a constant C such that for every configuration u and for q >0 

M^,,{T) > CT"-' . (10) 

Remark 3 It is important that the initial condition in Q is |0) and not an arbitrary state 
ip G £^(Z). In fact, could be an eigenstate of H^^ and hence not lead to any diffusion. 



In order to study the behavior of the Lyapunov exponent in the vicinity of the critical energy, 
that is, go beyond the trivial upper bound |7(-E'c + e)| < C let us define the transmission 
and reflection coefficients a^. and b"^ by 

MTi^+'M-^ V = a'^v + b'^Aj , ^ = ^ ^ • (H) 

Both are polynomials in e. As v is an eigenvector of all rotations, one has 

al = e^^'i , bl = . (12) 
Furthermore let us set e*''± = a!^/|a!^| so that rj^ = rj^-. Averages will always be denoted as 

(c±) = P+C+ +P-C_. 

Theorem 2 Suppose that (e^*''*) ^ 1 and (e^*''±) ^ 1. Then the Lyapunov exponent of a 
random polymer chain satisfies 



7(^c + e) - |i_(e^^.i)|2 + ■ (13) 
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Definition 2 A critical energy Ec G ^ of a polymer Hamiltonian is said to be of order r if 
\b'j.\ = 0{e^), but not = 0{e^^^) for both polymers. 

Remark 4 If is a critical energy of order r, tlien 7(-Ec + = C e^'' + 0{e'^^'~^^) for some 
non-negative constant C. Since 6!^ = 0{e), the order of every critical energy is as least 1, 
i.e. the Lyapunov exponent vanishes at least quadratically. Generically the order of a critical 
energy is 1 (this is the case in the dimer model). In the latter case, more explicit formulas for 
the coefficient in ()13|) invoking only the values of Tf" and OeT^" at the critical energy can easily 
be written out. A comparison of (fT^ with a random phase approximation is made in Section 
14.71 Finally, let us note that (fT^ also proves positivity of the Lyapunov exponent close to Ec 
whenever the numerator does not vanish. 

Remark 5 The conditions (e^*''*) ^ 1 and (e^*''*) 7^ 1 in Theorems |21 and 121 below are linked 
to anomalies studied in |CKj. For the dimer model, the condition (e^*''±) 7^ 1 is verified if and 
only if A 7^ l/"\/2- This particular value already appeared in ^BGj . 

Theorem 121 shows how the localization length diverges at the critical energy. The next result 
concerns the asymptotics of the integrated density of states M (denoted IDS, its definition 
is recalled in Section EiH below). Let Af± and A/!^ denote the absolutely continuous IDS and 
their densities associated to the models with L-|--periodic models composed of only one of the 
polymers. By definition of a critical energy, Af^{Ec) > 0. 

Theorem 3 Suppose that (e^"'*) 7^ 1. Then the IDS of a random polymer chain satisfies 

U(E.^,,.(l^^^,^l±mm^O(^). (14) 

The theorem states that M is linearly increasing at Ec so that there are many states in the 
vicinity of a critical energy. The spreading of these states is quantitatively nicely characterized 
by the diffusion exponents, namely the power law growth exponents of the moments M^^q(T) 
defined in (0) above: 



Here lim^ denote the superior and inferior limit respectively. The main result is the following: 
Theorem 4 Suppose that |(e^*^±)| < 1. Then -almost surely 

£ 1 - ^ • (16) 

It is interesting to compare Theorems IH and Q The latter implies for all configurations a 
weaker lower bound in (fTHI) of the form 1 — One can construct configurations uj with slower 
transport than in ()16|). Therefore - seemingly paradoxically - typical random configurations do 
not lead the slowest possible transport for this model. 
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Finally it is worth mentioning a large deviation result here. The IDS and Lyapunov exponent 
are both averaged quantities describing the behavior at the infinite volume limit. Given their 
asymptotics7Vr(Ee+e) = M{E^)+M\E^)t+0{t^) and7(Ee+e) = Ct^ + 0{t^) (here C = if the 
order of is p > 1), one therefore expects that typically (w.r.t. P) the following holds for the 
finite (but sufficiently large) size Hamiltonian H^^ jq found by restricting H^^ to ^^({0, . . . , N — 1}) 
(with Dirichlet boundary conditions): H^^^ has ciV^/^ equally spaced eigenstates in the interval 
[Ef. — A^~^/^, Ef. + A^~^/^] which are all spread out over the whole sample. Here we give an upper 
bound on the probability of the set of atypical configurations for which the average metal-like 
behavior of the eigenvalue spacing does not hold. This result shows on which scales there is 
strong level repulsion. 

Theorem 5 For every a > there exist c > and C < oo such that for all N E N there are 
sets Qn{(^) C Q satisfying 



P(n^(a)) = 0{e 



-cN° 



such that for every configuration uj in the complementary set VLiq{aY = Q\Qi\f{a) the following 
statement holds: the interval [Ec — N~^^'^~°', E^ + N~^^'^~"] contains of the order of N^^'^"" 
eigenvalues of H^^ j^i which are equally spaced and have eigenf unctions spread out over the whole 
sample, namely adjacent eigenvalues E and E' satisfy 

^ < \E - E'\ < ^ , (17) 



CN - ' ' - N 

and for all normalized eigenf unctions ip of H^^jq it holds that 

< mk-i)\' + \m\' < ^ (18) 

forO < k < N -1, where ^{-l) = ip{N) = 0. 

Outside of the interval [Ec — A^"^/^"", Ec + A^~^/^~"] we expect Poisson statistics. It seems 
unknown what the level statistics is like on the boundaries of this interval, but it is possibly not 
of the Wigner-type. 

Theorems 12] and El are proved through the perturbation analysis of polymer phase shifts and 
action multipliers (essentially, appropriately modified Priifer variables). The key for the proof 
of Theorems EJ and is Theorem which states that with high probability norms of the transfer 
matrices T^{k,m), 1 < m < k < N for energies in the interval [Ec - N'^^^-", Ec + iV-^/^-"] 
are uniformly bounded. This theorem is proved in Section El by establishing large-deviation 
estimates for random Weyl-type sums defined in terms of polymer phase shifts. 

Let us conclude with a brief remark about the one-dimensional Anderson model in the weak 
coupling limit, namely Hx ,^ = Ho + XV^j where Hq is a periodic operator, K; the usual Anderson 
potential and A a (small) couphng constant. Pastur and Figotin |PFj showed (in the case 
where Hq is the discrete Laplacian) that away from band-center and band edges of the periodic 
operator, the Lyapunov exponent grows quadratically in A. The large deviation results and 
dynamical lower bounds presented here transpose in order to show that almost surely 
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T>0 

Hence the presented techniques allow to study in a very detailed way the metal-insulator tran- 
sition driven by either the disorder strength or the sample size. This transition appears at a 
single energy, the critical energy, in the polymer models studied here. 



3 Brief review of basic formulas 
3.1 Transfer matrices 

Let (t{n))nez be a sequence of positive numbers and (f a sequence of real numbers. As 

in (j2I) they define a Jacobi matrix H acting on £^(Z). Given an initial angle 6*° G M and a 
complex energy 2; G C, let us construct the formal solution {u^{n))nez by 

-t{n + l)u''{n + l) + v{n)u''{n) -t{n)u^{n-l) = zu'{n) , (19) 
and the initial conditions 

t(0) u%0) \ _ f cos(e°) 
u'{-l) y V sin(^°) 

Using the definition (jS)) of the single site transfer matrices T^^t, the transfer matrix from site k 
to n is introduced by 

k 

T'{n,k) = Y[ T^(i)-zm ■ 

l=n--l 

It allows to rewrite the (formal) eigenfunction equation (fT^ as 

Note that the transfer matrices satisfy the transitivity relation T^{n,k) = T^{n,m)T^{m, k). 
A direct inductive argument then shows that, for ^ G C, 

r^+^(n, k) = r^{n, ^) - C E ^'^^(^' ^ + ^) ( J ) ^^(^' ^) • (21) 

Taking the norm of ()21|) . estimating the r.h.s. and consequently taking the supremum over 
< k < n < m leads to the following perturbative result which in a slightly different form is 
given in |Sim2j (see also |DTj .) 

Lemma 1 Suppose 

sup \\T'{n,k)\\ < C D = sup ^ . 

0<k<n<m 0<l<m-l \tl\ 
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Then, as long as CDKlm < 1, 

sup \\T^+<{n,k)\\ < ^ . 

0<k<n<m i — UUlC^lm 

3.2 Free Priifer variables 

Let now E & M. and be given by (fT^ . The free Priifer phases 6^'^{n) and amphtudes 
R^'^{n) > are now defined by 

^ ^""^ ^ sin(^^0'^(n)) ) ~ \ u''{n -I) ) ' ^ ' 

the above initial conditions as well as 

-1 < ^^o,E(^ + 1) _ ^o,£;(^) Stt 

Note that the ^^-dependence of the Priifer variables is suppressed. 
Lemma 2 

i?°'^(n)2 9z,e°'^H = <^ (23) 

I -T^ilnu'^ilf ^f n<0. 

Proof: From the recurrence relation (fTTHl and the definition of 6^'^{n) one gets 

cot(e°'^(n)) = -t2(n- l)tan(^°'^(n- 1)) +t;(n- 1) . 
Differentiation leads to 

= '^It^y^^'y^ d^e^'^'in - 1) + sin^(g°-^(n)) . 
cos^(t^'^'^(r2 — 1)) 

Multiplying with R^'^{nY and using the definition of R^'^{n) and 9^'^{n) gives 

/2°'^(n)29^j^°'^(n) = - - 1) + u'^iji - 1)^ . (24) 

The above deduction of has used that u^{n — 1) 7^ 0. If u^{n — 1) = 0, then one may 
deduce fj24|) in a similar way from 

mnPr ^^ COt(^°'^(n - 1)) 

tan(^°'-^(n)) - V V 77 



-t2(n - 1) + {v{n -l)-E) cot(0O'^(n - 1)) ' 
The lemma now follows by iterating (jSl)). □ 

Note in particular that (j23|) implies that dEO^'^in) is strictly positive for n > 2 and strictly 
negative for n < —2. Furthermore, it follows from elementary considerations for transfer matri- 
ces that there are constants Ci and C2 such that 

< Ci < |9£;^°'^(n)| < C2 < 00 , (25) 
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where Ci and C2 can be uniformly bounded away from and 00 as long as \n\ > 2 and 
the quantities \n\, E and vci8LX.\k\<\n\{\v{k)\,t{k)^ l/t{k)} remain bounded (where only the lower 
bound requires \n\ > 2). 

Let Iljv be the projection on ^^({0, . . . ,N — 1}) and denote the associated finite-size Jacobi 
matrix by = UnHUn. As has Dirichlet boundary conditions, let us choose u^{—l) = 
and t{0)u^{0) = 1 as initial conditions in the recurrence relation ()19j] . This corresponds to an 
initial Priifer phase 6^ = 0. The formal solution then gives an eigenvector (and E is an 
eigenvalue of Hn) if and only if t{N)u^{N) = R°'^{N) cos{e^'^{N)) = 0, that is ^"'^(A^) = | 
mod TT (note herefore that u^{0) 7^ for any eigenvector of Hn). 



One checks iteratively for all n > that 



[n 



> for sufficiently close to —00 and 



limE- 
limE- 



[n] 



l)/u (n) = 0. This and the definition of the Priifer phases implies that 
= for all n > 0, which one uses for n = N . As 6^'^{N) is monotone increasing 



in E, it follows that the jth eigenvalue Ej of (counted from below Ei < E2 < 
satisfies 

^o.^.(iV) = | + vr(j-l), 9' = 0. 
This oscillation theorem implies immediately: 



< E 



N, 



(26) 



- 9^''^{N) - # {negative eigenvalues of {Hn - E)} 

TT 



1 

< - . 

- 2 



(27) 



3.3 Modified Priifer variables 



Let us fix M G SL{2,R). Set eg 



cos(^) 
sin(^) 



. Define a smooth function m : 



with 



m{9 + tt) = m{9) + n and < Ci < m' < C2 < 00, by 



r{9)a 



m{e) 



Me, 



r{9) > 



m(0) G [-7r,n) 



X 



for initial condition 9 (0) 



Then the M-modified Priifer variables {R^ (n) , 9^ (n)) G 
9 = m{9^) are given by 

^^(n) = m(^°'^(n)) , 

and 

R^{n)cos{9^{n)) \ _ ( t{n)u^{n) 
R^{n) sm{9^{n)) ) ~ \u^{n-l) 

where the dependence on the the initial phase is again suppressed. Bounds of the form 
(j25j) also hold for an M-modified Priifer phase 9^{n) because 9^{n) = m{9^'^{n)) leads to 
(minm')|<9ij^°'-^| < |(9e^^| < (maxm')|9ij^°'^|. Furthermore, as \9'^{n) - 9^''^{n)\ < 27r, ^ 
implies that for the choice 6' = m(0) 



(28) 
(29) 



1 



TT 



9^{N) — 7^ {negative eigenvalues of {Hn — E)} 



< 



(30) 



The goal to have in mind when choosing M is to make the M-modified transfer matrices 
as simple as possible so that the M-modified Priifer variables are easy to calculate. Whenever 
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E is in the spectrum, the most simple matrix to obtain is a rotation. Anything close to it can 
then be treated by perturbation theory. This is the strategy followed for the random polymer 
model below where M is chosen as in ((7j) 

Example: Let us consider an L-periodic Jacobi matrix H. If G M is in the interior of the 
spectrum of H, there exists a matrix M (depending on E, of course) such that MT^{L)M~^ = 
where i?^ is the rotation by an angle t] = r]{E) obtained in accordance with the definition 
(f^H|) . The M-modified Priifer variables are then simply given by {R^{kL),6^{kL)) = {l,krj) 
and the IDS is Af{E) = r]{E)/{LTT). 

3.4 Covariant Jacobi matrices 

Let {Q, T, Z, P) be a compact space Q, endowed with a Z-action T and a T-invariant and ergodic 
probability measure P. For a function / G L^(f2,P), let us denote E(/(ci;)) = f dP{uj) f{uj). 
A strongly continuous family {H^)^(zfi of two-sided tridiagonal, self-adjoint matrices on £^(Z) is 
called covariant if the covariance relation UH^U* = H^uj holds where U is the translation on 
£^(Z). H^^ is characterized by two sequences {t^{n))nez and {v^{n))nez such that Q holds. 

The IDS at energy E eM. of the family {H^)^^q can P-almost surely be defined by |PFj 

^^iE) = lim ^ Tiixi-oo,E]{^NH^nr,)) , (31) 

while the Lyapunov exponent j{E) for G M is P-almost surely given by the formula 

7(i?) = lim llog(||rj^(iV,0)||) , 

where the transfer matrix Tj'^N, 0) from site to is defined as in Section I^TTl Both the IDS 
and the Lyapunov exponent are self-averaging quantities, notably an average over P may be 
introduced before taking the limit without changing the result |PFj . 

For each if^ let {R^ (n), 0^{n)) denote the associated M-modified Priifer variables with some 
initial condition, then according to (jHUj) 

Af{E) = hm i 1 E {e^{N)) . (32) 

While it is readily seen that 'y{E) > lim^_^^ j^Eilog^R^ (N)) , one may in general not get 
equality here as demonstrated by a counterexample in Section H?Tl This is due to the dependence 
of R^ (N) on the initial phase 6. The next lemma solves this problem by (continuously) averaging 
over 9. 

Lemma 3 For G M and any continuous [i.e. non-atomic) measure v on RP(1) = [0, vr) 

7(i?) = jhn^ ^ 1 dv{Q) E (log(i?f (AT))) . (33) 
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Proof: As ||7^^(A^, 0)em,-i{e) || = ||^em-i(6»)|| ||^ ^^eS{N)\\R^i^)y ^ change of variables and 
elementary estimates show that it is sufficient to show that 7(-E) is equal to 

i^ooJl J M^)E(log(||rj^(iV,0)e,||)) 

for any continuous probability measure u. This is easy to see if 'y{E) = 0, thus we now assume 
that 'j{E) > 0. Suppose the contrary, that is there exists a v such that 

Jhn I du{e)E{\og{\\T^^{N,0)ee\\)) < i{E) . 

By Fatou's lemma this implies that J diy{e)'E{\imN^^ j^\og{\\Xf{N,0)eg\\)) < 7(E). Be- 
cause for a.e. a; the limit inside of the expectation is equal to either 'y{E) or — 7(-E) by 
Oseledec's Theorem, there has to exist a set £^ C [0, vr) x of positive u ® P-measure such 
that limAT^oo ;^ log(||'^(^) 0)e6)||) is equal to — 7(-E) for all {0,lj) G £. Hence there exists an u 
such that the set {9 G [0, vr) | E eigenvalue of H^{9)} has positive u-measure where H^{6) is the 
half-line operator with ^-boundary condition. As v is continuous, this set has to contain at least 
two distinct points. This is in contradiction to the fact that the difference equation H^u = Eu 
has, up to constant multiples, at most one square-summable solution at +00. □ 



4 Asymptotics of IDS and Lyapunov exponent 

Generalizing the strategy suggested by Pastur and Figotin |PFj . this chapter is devoted to the 
calculation of the asymptotics for the IDS and the Lyapunov exponent near the critical energy 
of a random polymer model, that is the proof of Theorems El and El The techniques of jC^Sj 
allow to treat also the case of strongly mixing (instead of random) configurations of polymer 
chains giving similar formulas, containing a correction factor given by the Fourier transform of 
the correlation function. No further details are given here concerning this generalization. 



4.1 Random polymer chains 

For sake of completeness, let us briefiy indicate how to construct (fi, T, Z, P) for the random 
polymer Hamiltonians defined in Chapter |21 Let f2o be the Tychonov space of two-sided se- 
quences of signs. Set VL± = {uo E VLq\ujq = ±} x {0, . . . , L± — 1} and = U Now 
T : Q ^ Q is defined by 

r {uj,i + i) if / < - 1 , 

T{u;,l) = \ 

[ (To^,0) if / = L^„-1 , 

where Tq is the left shift on Qq. Now for any set A± C Qq of codes all having ujo = ±, one sets 
for all / G {0 . . . L± - 1} 

P{{{u,l)eQ\u}eA±}) =- -^^^^^^ 
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where Pq is the Bernoulh measure on f2o- It can then be verified that P is invariant and ergodic 
(the latter by mimicking the proof for (f2o, To, Pq))- Random hopping terms and potential 
are then given by = (. . . , ta;o, 4i, • • •) and v^^^j) = (. . . , -O^g, -O^^, . . .) with choice of origin 
t{uj,i){0) = 4o(0 and t;(a;,o(0) = v^oiO- This leads to the covariant family (-ff(a;,o)(<^,OGf7 of 
Jacobi matrices. It is this family which is refered to as (H^^) in Section |21 and, in particular, in 
Theorems |21 to El 

According to Section 13.41 the Lyapunov exponent satisfies 

liE) = jhn^^E (log(||7;f,)(iV,0)||)) = Jim 1 log(||7^f,)(iV, 0)||) , (34) 

for P-a.e. {u, I) G f2. Here E = J dP. On the other hand, there is also a Lyapunov exponent 
associated with random products of the unimodular matrices Tj': 

7o{E) = hm I Eo (log (||T5(A;, 0)||)) = lim i log (||Tf (A:, 0) ||) , (35) 

K- — ^■OO fx, fc~^00 fx, 

for Po-a.e. u E Qq and Eq = / c/Pq. To compare 7(-E) and 7o(-E), let CIq be the full measure set of 
those uj E flo such that holds and also X]z=o^ ^i^J^ ~^ (-^±) as A; — * oo. For u E Qo^i is easily 
seen that lim^v^oo ^ log(||7^Jo)(^. 0)11) = lo{E)/{L±). Since P{(cu,0) \ to e n^} = > 0, 

one concludes from ()34|) that 

liE) = loiE) . (36) 

While 7o is not defined through a covariant operator, it follows by the same argument as in 
Lemma 121 that for any continuous measure u on [0, vr) 

7o(i?) = lim I [ du{9)Eolog{\\MT^{k,0)M-'ee\\) . (37) 

Counterexample: The continuity condition on u in Lemma 01 cannot be weakened as shows 
the following example. Consider the polymer model with L± = 3, t{l) = 1 for all / G Z, and 
v+ = (|,2,0) and V- = (— 2, 0), and choose M = 1. For _E = it is easily seen that 
T^Gt,/2 = =F|e7r/2 and thus ^ log ||T°(A;, 0)e^/2|| = — | for all u and k, while 7o(0) = |. Hence a 
measure having an atom at 6' = ^ will not satisfy (jTfj) . This also provides a counterexample to 
Lemma El with (f2,T, Z, P) as above. For this one uses that the event {^^[^(^(O) = 'Oa;o(0)} has 
probability 1/3 in Q. 

4.2 Polymer phase shifts 

For M given by ((71), let the polymer action multipliers p'^^{9) and the polymer phase shifts S^_±{9) 
be the M-modified Priifer amplitude and phase for the L-t-periodic polymers with initial phase 
6^ at and evaluated at L± (i.e. over a single polymer (t_|_,-0_|_) and (£_,{)_), respectively). By 
definition of the modified Priifer variables, this means 

P±{0)es,Me) = MTt+'M-'ee (38) 
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for all G M. The iterated polymer phase shifts are then denoted by 



From ((Tj) it follows that (independent of 9) p\{9) = 1 and ?7± = 50,^(6*) — 9, at least up to a 
multiple of 27r which is hereby fixed. The former readily implies that 7(-E'c) = 0. To study the 
Lyapunov exponent in a vicinity of E^, iterate in order to deduce 



7V-1 



log(||MTj=+^(iV,0)M-^e,||) = log (p^^^^J^))) , (39) 
which combined with fj36|) and (jHTf) gives 



1 1 ^'^ r 

7(i?e + 6) = -—Jim du{9)Bo{\ogip^SU9)))) . (40) 

To also express the IDS in terms of the polymer phase shifts, let {n(^i^^i)^k)kez be the sequence of 
lower polymer nodes for a given (cu, /) G fl, i.e. the integers determined by V(^i^^i){n(^i^^i)^k) = v^^^^O), 
for any choice of v. For G N, let ?T-(a;,«),A: be the polymer node closest to A^. Since S^^{9) — 9 
is a rotation number for a matrix which arises from H^i^^i^^n by a perturbation of rank bounded 
by Cmax{L_,L+}, it follows that \9f^^^\N) - {Sl^{9) - 9)\ < Cmax{L_,L+} uniformly in 

9. Thus it follows from (jH^ that Af{Ec + e) = limjy^oo ~j^{<Se.uji.^) ~ ^) almost surely and in 
expectation. Since k/N ^1/ {L±) almost surely as — > oo, this implies that 

M{E, + e) = lim Ieo{S^J9)-9) 



^ 1=0 



(41) 



4.3 Calculation of phase shifts and action multipliers 



The aim of this paragraph is to calculate the polymer phase shifts and action multipliers needed 
in (PT|) and (jin|) in terms of the transmission and reflection coefficients defined in fjllj) . Because 
det(MTf^+'M-^) = 1, these coefficients satisfy 

= 1 . 

A further short calculation shows that 

p^(^)2 = l + 23?e(a^6^e2*^) +2|6y% (42) 

and 

g*(5,,±(e)-e) ^ a± + b^e "^'^ 
\a% + b%e-^'^\ 
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Now using the phase of a!^, 

a'^ = e^^i + 0(|6^n . (43) 
This leads to the following expansions: 

log(p^(^)2) = 2 3fJe {a^ble^''^) + \b'J^ - {{a'^b'^fe^'^) + O {{b^J^) , (44) 

and 

^Hs.Me)-e) ^ e^^^i + fo^e^^i e-^*^ - fe^e^^^i e^^^ + O . (45) 

4.4 Oscillatory sums 

Proposition 2 Let c± G C, j = 1, 2, and set 

N-l 

1^0, e) = Eo(/i,^(e,6)) , Ii^^{9,e) = c., e^^^'^-^^) . 

/=o 

Let e be sufficiently small. If (e^*^''*) ^ 1, then lij{9,e) = 0{N\b'J, 1). // (e^*^''*) 1 /or &ot/i 
J = 1,2, 

lU0,e) = N{c^) ^_^^2rni^ +0{N\bl\', 1) . 
Proof: Since S'^lJi^O) = Se^ujiiSj, ^{6)) and 5' ^(6*) is independent of ui, one gets 

N-l 

11,(9, e) = (e^^^i) J^_,(e, e) + (c±) e'^' + (c±) Eo (e^^^-^^^-^'^)) - e^<''".+'5-(''») . (46) 

1=1 

Equation (gSD shows that e^*'^.,.,^ _ g2*(^^^+e) _ od^^j). As /^(^,e) = 4_i(^,e) + 0(1) and 
(e^*''±) 7^ 1 by hypothesis, one can solve for l]q{9,e) which directly implies that llf{6,e) = 
0{N\b''^\ 1). Along the same lines, l'^{9,e) = 0{N\b'±.\, 1). Now insert the expansion in 
Due to the above, the oscillatory terms in that formula are then of order 0{N\b'^_^\'^). Thus 
only the non-oscillatory term on the r.h.s. of (jl^ gives a contribution to the leading order. □ 

4.5 Asymptotics of the IDS 

Proof of Theorem H Formula (03) leads to 5,,±(^) - 6 = r]± + ed± - e'^m{c±e^'^) + 0(6^), 
where d± = (9e?7^)|^_g and c± = id^b'l^)\ e^^^ . Inserting this in yields 

X{E, + e) = ({v±) + e - e hm i Eo f V c^,e'<-^A + Oie' 

'7l[L±) \ fc^oo fc > ^ — ^ ' 



. «=0 
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By Proposition 121 the expectation of the oscillatory sum is of order 0{k\b'^_^\, 1) and thus 

^fiEc + e) = ((r^±) + e(rf±) + 0{e')) . (47) 

Setting p+ = 1 and p+ = yields that in particular J\f±{Ec + e) = ^;r^(r7± + ed± + 0{e^)), 
allowing to identify J\f±{Ec) = Tj±/TiL± and A/'4(-Ec) = d^/nL^. Using this to insert for rj± and 
d± in (jTTj) completes the proof. □ 



4.6 Asymptotics of the Lyapunov exponent 

Proof of Theorem 121 Replacing ()44j) into ()40|) shows that 2{L±) 'j{Ec + e) is, up to corrections 
of order equal to the i/- average of 



( N-1 \ / N-1 

1=0 / \ 1=0 

By Proposition 121 the first oscillatory sum has a contribution of the order (which is 

given there) while the second oscillatory sum is of order (9(|6^|^) and can hence be neglected. 
Therefore one obtains 



7(^c + e) 



{L±) 



\bl\^)+^e 



1 - (e^^^i) 



(48) 



It can be directly verified that the given leading order term vanishes if either p+ = or p+ = 1, 
which also follows from the fact that in this case H^^ is a periodic Jacobi matrix, whose Lyapunov 
exponent vanishes in the interior of its spectral bands. Next rewrite ()48|) as a fraction with 
common denominator |1 — (e^*''±)p. Since p_ = 1 — p+, the numerator is a polynomial of degree 
at most 3 in p+ vanishing at p+ = and p+ = 1. Elementary but lengthy algebra shows that 
moreover its third derivative vanishes identically. Calculating the first order derivative allows 
to conclude. □ 



4.7 Comments 

A random phase approximation consists in supposing that the incoming phases Sl^{9) in 
each summand of (pn|) and (^1]) is completely random, that is distributed according to the 
Lebesgue measure. It can easily be checked that one actually obtains the correct answers for 
the derivatives of both the IDS and the Lyapunov exponent at the critical energy within this 
approximation. However, the lowest order non-vanishing term in the Lyapunov exponent is 
0(|6!tP) and the random phase approximation gives together with the expansion ()44|) that 
^{Ec + e) ~ (|6!^p)/(2 {L±)), namely only the first term in (PH|) . As we shall argue now, the 
second contribution is due to the presence of correlations (or memory) in the family of discrete 
time random dynamical systems (5^,-1-, ]RP(1), f2, P)^^^. 
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It is a result of Furstenberg |Furj (his hypothesis can be checked here) that for each e 7^ 
(small enough) there exists a unique invariant measure on MP(1) satisfying 

j duM f{6) = I du,{e) (/(5,,±(^))) , / G C{RP{1)) . 

For e = 0, one invariant measure is given by the Lebesgue measure (it is unique if 77+ — rj^ is 
irrational). For finite e, iteration of the invariance property and the Proposition |21 implies 

Jdu.{e)e-^ = |c^.,(^)E(^l^'e-<^(^)j = ^^^^ + ^d^^^) • 

Similarly J dv^{9) e^^^ = 0{\b'^\). These facts express the deviations of the invariant measure 
from the Lebesgue measure and hence from the random phase approximation. Moreover, for 
small enough e 7^ the invariant measure z/^ is known to be Holder continuous |BH p. 161] so 
that one can use it in (j4U|) . Hence 

7{Ec + e) = / d^O) {\og{pWm ■ 

Developing log(p!^(^^)^) as in ()44|) then also leads an alternative proof of ()48|) and the second 
contribution in PHj) is indeed due to the correlations as claimed above. Finally let us point out 
that higher order terms in e can readily be calculated, under adequate (weak) hypothesis. 



5 Large deviation estimates 

Using elementary estimates on the boundary terms M and M^^ in ()39|). as well as (fT^ . and 
the expansions ()43|) and ()44|) . one obtains that for all < m < /c < iV 

k-l 

logf||Tj=+'^(A;,m)||') =26 sup J" c^^ e^^-'-'^ + 0{N6\l) , (49) 
^ ^ em.) 

where c± = e^^^{dsb^^)\s=o. If the order of the critical energy is 1, then c± = 0{1). In order to 
prove the delocalization results, it is necessary to show that the l.h.s. of (jl^ is of order 1 as long 
as OlNS"^) = 1. Therefore one needs to show that sums like I^js{d,S) defined in Proposition |21 
are with high probability of order y/N for 0{N6'^) = 1 and \k\ < N. These random Weyl sums 
can be thought of as a discrete time (variable N) correlated random walk in the complex plain, 
the correlation being due to the presence of the dynamics Ss^±. For the present purposes, it is 
sufficient to show that this sum actually behaves as a random walk on adequate time scales. 
Hence let us introduce, for every 6, 6, 

n%{a,6,e) = {cj G 1^0 I 3 A; < such that |/^^fc(^,(5)| > A^"+^/2} . (50) 
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Theorem 6 // |(e^*''*)| < 1 and a > 0, there exist constants Ci and C2 such that for all 6, N 
and 6 with N6^ < 1: 

PoinUa, 6,6)) < Cie-^^^" . (51) 

The proof of this estimate wiU be given in Section 15.11 First, let us deduce the following 
consequence: 

Theorem 7 Let |(e^*''='=)| < 1 and a > 0. Then there are c,c' > 0, C < 00 such that for every 
N eN, there exists a set Qi^{a) C Q satisfying 



P{nN{a)) = 0{e 



-cN° 



and such that for every configuration {uj,l) in the complementary set Q^la)'^ = Q\Qn{c(), one 
has 



m] 



< C , 



for alio <m<k<N and all \6\ < A^-"-i/2^ |k| < c'/N. 

Proof: In order to estimate the norms of the transfer matrices using the Weyl sums, note that 
for any 2x2 matrix A, 

\\A\\ = sup ll^eell < \/2 max II Aeell . (52) 

Set n%{a,S) = Q%{a,6,0) UQ%{a,6,^). Then, combining Theorem with ^ and ^ 
as well as the fact that T{k,m) = T{k,0)T{m,0)~^, one deduces that for all u G r2^(a,5)'^ 
with \5\ < A^~"~i/2, norms of the transfer-matrices T^''~^^{k,m) are uniformly bounded by a 
constant, not dependent on 6. 

Now let QN{a, 6) = {{u},l) e Q \ uj e ^l%{a, 6)}. It follows that 

P{n^ia,S)) < + Po{n%{a,S)) < Cse-^^^". 

Elementary estimates (based on uniform bounds on norms of transfer matrices over blocks 
of length no more than maxL±) imply that 

||7^5f(fc,m)|| < C", 

for all (to, I) e fi7v(a, 5)^, \S\ < A^-"-i/2 and < m < A; < (in fact this holds for m, k up to 
the A^-th polymer node). Set e = A^^"^^/^. The theorem then follows from Lemma ^ by taking 
^^(«) = UM^^(«.^e/Ar). □ 
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5.1 Correlation bounds: Proof of Theorem [SI 

Lemma 4 Let n = |(e^*''±)| < 1. Then there exists a centered complex random variable X{uj) 
depending on ui, . . . ,ujr such that 

e'^^l-^'^ = Xiuj)e''' + Oir6,K'') . 
Moreover, \X{uj)\ is uniformly hounded by 2. 

Proof: Let us set Kr{uj) = exp(2z ^^^-^ 77^^). Note that |Eo(Kr(ci;))| = k''. Iteration of 
^2iSs,±(0) — g2i(r?±+e) _j_ Q(^§^ centering the random variable Kr{uj) shows 

g2.5J,„(e) ^ ^^^^^^ _ (^^2^v±yry^e ^ ^(^r^ ^ 

as claimed. □ 

Proof of TheoremEl Let r be the smallest integer larger than log(A^^"^^/^)/ log^n). Applying 
Lemma m to each term (except the first r terms) of the sum I^ i:{9, S) shows that 

fc-i 

= Y. ^-^^r ^) e'^^^-^'^ + 0{rk5, kK\ r) . (53) 
1=0 

(Here the identity Sl'\J{9) = <Ss^Tlui'^s,uji^)) used. Recall moreover that TqLU is the /-fold 
shift of u.) Under the hypothesis of the theorem and because of the choice of r, the error term 
in ()53|1 is 0{N^^'^ log N). Thus it is sufficient to prove probabilistic estimates of the appearing 
sum, which will be denoted by Zk{u!, 6). 

In order to decouple the correlations, divide {0, . . . , A; — 1} in 2R pieces Iq, ■ ■ ■ , hn-i of equal 
length [k''], where i? = [k/ {2[k''])], i.e. J, = {s[k% {s + l)[k''] - 1}, s = 0, . . . ,2R - 1. Here [x] 
denotes the largest integer smaller or equal to x. This excludes ck°' terms which in the following 
can be absorbed in the error. Set for j = 0, 1: 

R-l 

s=Q leis 

Thus 

Zkiuj, 6) = Zl{uo, 6) + Zl{uo, 6) + 0{k^) . (54) 

The random variable Ys{ijj,5) satisfies uniformly |ys(co',5)| < Ci fc". If denotes the averaging 
procedure (conditional expectation) over all random variables uji for / > s, then LemmaEJimplies 

E,[fe«]+i(F,(^,5)) = 0. 

In the following estimates, real and imaginary parts of Z-j^{u,6) are treated separately, but 
in exactly the same way; hence one may suppose that Z-'j^{uj,5) and all the summands therein 
are real. For A > and /? > 0, the Tchebychev and Cauchy-Schwarz inequalities imply 

Po({^Gi7o|^fc(^,5)>A}) < e-^^ Eo(e^^^(-'^)) < 6"^^+^'^^" max Eo(e2^^«(-'^)) . 

i=o,i 
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Now if -1 < F < 1, by convexity 2 6^^^ < (1 - Y)e-^ + (1 + ^)e^. Thus if F is a real 
centered random variable, 

E(e'^^(-)) < (e-^ + e^)/2 < e""'^ . (55) 

One may assume that [A;"] > r — 1 and k > (otherwise it is trivially true that /^{O, 6)\ < 

Thus Zj^_^{uj,S) does not depend on the cui with / > {2{R - 1) + j)[A;°] - 1 and a 
rescaled version of can be iteratively applied to the conditional expectations, leading to 



Eo(e2^^«(-'^)) < Eo (E(2(^-i)+,)[,.]_i(e^^^^(--)+^(--^))e^^^«-(-'^)) 

<- g(2cifc"/3)2/2g^ J^g2/3Z^_^(a;,<5) j 

Choosing /5 = A/(2c2A;^"''") and proceeding similarly for {u E QqI Zk{uj,6) < —A} thus shows 
(after recombining real and imaginary parts) 

Po({.. G Qo I 1^.(^,5)1 > A}) < 4e-^^/(^-^^"")+^ . 

Using this estimate for A = A^"+^/^ and renormalizing the constants in order to compensate for 
the error terms in ()5Hj) as well as for summation over k concludes the proof. □ 



5.2 Eigenvalue distribution in the metalic phase 



Proof of Theorem IHl Let us fix a configuration (cu,/) G fiAr(a) (see Theorem [7j) and suppress 
its index. Using pn|) and the fact that the norm of a transfer matrix is equal to the norm of its 
inverse yields ||T^'=+''(fc, 0)||-^ < R^'^''+^{k) < ||T^'^+'^(A;, 0)||. Theorem [Z| therefore guarantees 
the existence of a constant C such that for < < and for — A^^^^^/^ < ^ < jv^°^^/^, 

1 < RO,E^+s^f^Y < 1 < |n^=+^(A;)|2 + |M^^+'^(fc-l)|2 < C . 

This readily yields (fTH|). Now by Lemma |2l 

^ < ai5 0°'^^+^(iV) < NC . 

Upon integration, one deduces 

^\E-E'\ < |5°'^(iV) -^°'^'(Ar)| < NC\E-E'\, 

for all E,E' G [E^ - iY^^-^/^ + N-'^-^/\ The oscillation theorem discussed in Section 
now gives (fTTj) . □ 



20 



6 Lower bound on dynamics 



The deterministic part of the argument presented in this section follows jPTj . Let us return to 
the simplified notation from Section El and write uo instead of (cu, /) since based on the results of 
SectionElthe value of / will not influence the considerations. Let us begin with some preliminaries 
and introduce the Green's function 

1 



Gl{n) = i n 

Note that 



-t^{n + 1) Gl{n + 1) + {vM -E,-z) G'M - t^{n) G'^{n - 1) = 6n,o • (56) 
Using transfer matrices, one now has for n < 0, 

G^(n-l) J ~ ^-^^'U^!^ G'^(-l) J ' ^ ^ 



while for n > 1, 



Un)GM \ _ ^^ f ^..(1)^^(1) \ ..ox 

G^in-l) ) - ^-^^'1^^ Gm ) ■ 
The following identity is well-known: 

M.,,(T) = i 1 ^ / rfE \GUn)\' , z = E + ^. (59) 

Proof of Theorem |3] For given a > let c, c' > and C < oo be the constants form 
Theorem and choose N = [c'T] and e = A^"^/^"". By Theorem [3 there exists ^^^(q;) C f2 
with P(fi^(a)) = 0{e-'^") and such that for uo G fi7v(a)' one has \\%f'^^^'^^{n, 1)|| <G for 
all \5\ < A^-"-i/2 and n<N. 

For such u, because of the uniform bounds on the matrix elements t^{n) and Vt^{n), for 
n = one of the three terms on the l.h.s. of has to be large. Suppose first that |G^(0)p + 
1^^(1)1^ > Cg > 0, then it follows from ^ and ||(TJ(n, = \\T^{n, 1)|| that 

max{|G:(n)|MG:(n-l)r} > 



According to the above, as long as 5 G [— e, e] the transfer matrices are bounded from above 
by C as long as n < [ci T], in which case at least every second |G'5(n)p is bigger than Gt/C"^. 
Replacing this into (|^. 



0<n<[ciT] 



for some constant Cg > 0. If, on the other hand |Gf,(— 1)^ > Cg > 0, then one gets this 
estimate in the same way, but based on (jFf|l instead of (fSH|) . This uses the fact that the analysis 
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of Section El can also be carried out for T^"^^ {n^Q) with negative n. A Borel-Cantelli lemma 
shows that a.s. > 1 — {\ — a) / q. Since a > is arbitrary, this finishes the proof. □ 

Proof of Theorem Q Follow the above argument by using the deterministic Proposition ^ □ 
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